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Abstract: Members of the transforming growth factor-beta superfamily, growth differentiation factor 9 (GDF9) and 
bone morphogenetic protein 15 (BMP15), have crucial roles in fecundity of sheep. Our previous investigation confirmed 
that the fecundity mutations of sheep presented in highly prolific White goat individuals of Guizhou province. To 
illuminate other polymorphisms in BmpJ5 and Gdf9 genes and the relationship of these mutations with function, we 
cloned and characterized the coding region of Bmp15 and Gdf9. Molecular models of BMP15 and GDF9 mature peptide 
of White goat were constructed based on the homology of human BMP7 experimental tertiary structure. Two exons 
encoded prepropeptide of 394 amino acids in BMP15 and 453 residues in GDF9, respectively. Apart from the FecX?® 
mutation (S991) in BMP15 and V79I mutation in GDF9 confirmed in White goat previously, other seven and three 
polymorphism sites were detected from BMP15 and GDF9 mature peptides, respectively. S32G, N66H, S99I/P99I and 
G107R in BMP15 could be important for the binding of dimer to receptors. Changes of P78Q and V79I in GDF9 might 
affect the binding of dimer to receptor type 1. Comparing the length of BMP15 and GDF9 prepropeptide in vertebrates, an 
increase in length of BMP15 presented along with the protein evolution from fish to mammal and the divergence of the 
N-terminus residues in matured BMP15 peptide might contribute to the sensitive control on the fertility of animal species 
with low ovulation rate. These findings gave a valuable explanation for the correlation of mutations in Bmp/5 and Gdf9 
genes with the control on fecundity of White goat and supported the notion that they were the pivotal factors in female 
fertility of White goat in Guizhou province. 
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HE: BMPIS#IGDFOHILEKA LB CTGFB) HAIR Rai, SARE ERA Bera VER, MAP 
KIMI Pei RAE A Ehe TRON A I. FEIT, VE IN LB SE 
ARNE Ro YT EDF AR Bmp 15M Gdf9 EARE PEK ZI KAR, HDI A LUE Bmp 1 SA GadfO Fe Al hi IX, 
BEAT THEME, DL ABMP7 F HAA aE A RARE SN A LAE BMP15A#IGDE9 FARK) = AER. EN ALLE Bmpls 
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TAIZ AERA. FETA, BMPISAKFANKAYS32G. N66H. S99T/P99TAIG1O7R RA HY RESY Ui) — BR RG RE ANI G 
GDF9 SATA ASP 78QAIV 79 L532 Mid — 5 Ps LS PRISE A, REESE RATE. X Bmp 15 Al Gdf Ala 0S H 
AB AL BI ASF ETT RKN, SER REAKIRI SEE LE A, BMP 15 HH SK eT I LR, LABMP15 
RFA ENE KENIA E. PRR N AE EBMP 15 Xt (ECE SA FL eT REAR TCI GR A 
oN ALE Bmp 15 FG fo ENR ot II ABE TAER, JPRS SAFE IN ALE ee PEK 
Bed TAT AS o 
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Gene insertions and deletions have been widely 
acknowledged to play an essential role in the evolution 
of genes in vertebrates (Gemmell & Slate, 2006). Most 
mutations in protein are associated with disease, only a 
small number of examples show the heterozygote 
advantage. The good example is Bmp15 and Gdf9 genes, 
which are identified to be related with female fecundity 
in several animals. Bmp15 and Gdf9 genes encode for 
bone morphogenic protein 15 (BMP15, also known as 
GDF9B) and growth differentiation factor 9 (GDF9), 
respectively. Both of them are members of the 
transforming growth factor B (TGFB) superfamily, in 
which over 35 proteins have been found with pivotal role 
in cellular growth and differentiation (Juengel et al, 
2004). These two oocyte-derived molecules, BMP15 and 
GDF, are closely related in their primary structures and 
share a nearly identical spatiotemporal expression pattern 
in the ovary of animals (Juengel et al, 2004). However, 
these two molecules are encoded by distinct genes and 
BmpI5 gene maps to chromosome X but Gdf9 to 
autosome 5 in sheep (Davis, 2005). Active peptides of 
BMPI15 and GDF9 contain 125 and 135 amino acids 
respectively (Davis, 2005; Galloway et al, 2000). 
Currently, it is hypothesized that BMP15 and GDF9 may 
form non-covalent heterodimers or homodimers to 
regulate the fertility of small ruminants (Hanrahan et al, 
2004). The dimers combine and activate two types of 
membrane-bound serine/threonine kinases on the 
membrane of target cells in ovary. BMP15 is reported to 
bind to ALK-6 and BMPR-II (Moore et al, 2003), while 
GDF9 integrates to ALK-5 (Mazerbourg et al, 2004) and 
BMPR-II (Vitt et al, 2002). Mutations in Bmp15 and 
Gdf9 genes have been shown to increase ovulation rate 
and ultimately litter size in sheep (Hanrahan et al, 2004). 
Five separate mutations in the mature Bmp/5 coding 
2005) with 
heterozygous advantage in sheep. One copy of the 
Inverdale (FecX’), Hanna (FecX"), Belclare (FecX®), 
Galway (FecX®) or Lacaune (FecX") allele increases 


region have been described (Davis, 


about 0.6 lambs per ewe compared to the wildtype 
(Hanrahan et al, 2004; Bodin et al, 2007). Conversely, 
homozygous carriers of Bmp15 mutation have small 
undeveloped ovaries and are infertile. Similar to BmpJ5, 
invalid Gdf9 gene also blocks follicular growth in 
homozygotes, resulting in sterility in sheep, and one 
mutated FecG” allele lends to increase ovulation rates in 
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heterozygous ewes (Hanrahan et al, 2004). Furthermore, 
ewes heterozygous for mutations in both Gdf9 and 
Bmp15 are fertile and the effects of these mutations on 
ovulation rate are additive. It appears that both of Bmp/5 
and Gdf9 are important in the control of ovulation rate in 
sheep. However, the relative importance of these growth 
factors differs in mice. Gdf9 is essential for normal 
follicular development whereas BmpI5 is not because 
Bmp15 knock-out mice are fertile only with the ovulation 
and fertilization of oocytes impaired. It suggests that the 
regulation of BmpI5 and Gdf9 genes on fertility is 
species-specific in mammals. In nonmammalian 
vertebrates, both of Bmp15 and Gdf9 have been detected 
in ovaries of hen (Elis et al, 2007; Johnson et al, 2005) 
and fish (Clelland et al, 2006; Liu & Ge, 2007) but the 
correlation with their fertility remains to be elucidated. 
Goat normally has a low ovulation rate, which is 
similar to sheep. The mRNAs and proteins of Gdf9, 
Bmp15 and BMP receptors have been detected in goat 
ovarian follicles at all stages (Silva et al, 2004). 
Polymorphisms of Bmp/5 and Gdf9 of White goat in 
Guizhou province had been confirmed in our previous 
investigations (Du et al, 2008; Lin et al, 2007). In 
addition, it is characterized by litter size with extreme 
variation among individuals from one to five lambs, 
which is consistent with the phenotype of segregating 
major genes. Therefore, White goat is a desirable model 
to study the diversity of these genes for control of 
Thus, to 
polymorphisms in Bmp/5 and Gdf9 genes and the 


ovarian _ folliculogenesis. illuminate 
relationship of these mutations with function, the coding 
regions of Bmp15 and Gdf9 genes of White goat were 
cloned and a comparison with other animals was also 


undertaken in the present study. 


1 Material and Methods 


1.1 Animal and samples 

Five prolific White goats chosen for this experiment 
were healthy and recorded to give triplets birth. They 
kindly provided by Qinglong goat breeding farm in 
Guizhou province of China. Genomic DNA was 
prepared respectively from the blood sample of goats as 
the protocol of E.Z.N.A. TM SE Blood DNA Kit 
(Omega) described. Electrophoresis on 0.7% agarose gel 
showed this material to be of high molecular weight. 
1.2 Amplification and cloning of Bmp15 and Gdf9 
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genes 

The sense and antisense primers were based on the 
conserved sequences of Bmp15 and Gdf9 genes in goat 
and sheep deposited in GenBank (Tab. 1). Genes were 
amplified by polymerase chain reaction (PCR) with 
primers (50 pmol each), 100 ng of genomic DNA, 2.5 
units of Ex Tag DNA polymerase (Takara), 0.2 mmol/L 
of each dNTP, and 10 uL of 10xreaction buffer adjusted 
to a final reaction volume of 100 uL with 3dH,O. An 
initial long (5 min) denaturation step at 94°C was 
followed by 30 s at 50—56°C and 2min at 72°C, and 
then 30 cycles of 1 min at 94°C, 30s at 52— 60°C, and 1 
min at 72°C. The reaction was stopped after a final 
extension time of 10min at 72°C. The size and purity of 
the PCR products were estimated by subjecting the 
samples to 1% agarose gel electrophoresis. 


PCR products were cloned into the pMD18-T 
simple vector according to the procedures of the PCR 
cloning kit (Takara) and transformed into the competent 
Escherichia coli cells TG1 {supE A(hsdM-mcrB)5 thi A 
(lac-proAB) F’[tra A36proA+B+laclqgZDM15]}. The 
positive clones were detected by colony PCR. The 
plasmid DNA carrying the PCR products was prepared 
by the procedure of E.Z.N.A. Plasmid Miniprep kit 
(Omega) and subjected to sequencing by Applied 
3730 automatic DNA sequencer. The 
sequencing was carried out by Shanghai Invitrogen DNA 


Biosystems 


Sequencing Service. To correct for potential errors 
arising from PCR amplification, sequencing was carried 
out on at least two clones, derived from independent 
PCRs. 

1.3 Sequence analysis and tertiary model 


Tab. 1 Primers for amplification of Bmp15 and Gdf9 genes in White goat of Guizhou 





Primer name Primers sequence (5'—3') PCR products (bp) 

Bmp15 exon] sense GATGCAAAGAGGACAATTTAGAAGACC 574 
antisense CCCACCAGAACAATATAGTATGATAACTC 

Bmp15 exon2 sense TGCAGGCTCCTGGCACATACAGAC 862 
antisense TCACCTGCATGTGCAGGACTGGG 

Gdf9 exon1 sense GGAAGAAGACTGGTATGGGGAAATG 461 
antisense CTGCTCCTACACACCTGCCGC 

Gdf9 exon2 sense TCTTCTCAGGAACCTTTCCATCAGT 995 


antisense 


Searches for homologous DNA and protein 
sequences were conducted with the BLAST software 
against the nonredundant GenBank database (http:// 
www.ncbi.nlm.nih.gov/blast/blastn/). Phylogenetic an- 
alysis was performed with MEGA4.0 program using 
Neighbor-Joining method. The confidence level for each 
branch was estimated by interior branch testing with 
1000 random replications using Clustal W method. 
Predictions of the structure were performed online with 
the Phyre program (Bennett-Lovsey et al, 2008) and 
deepview/Swiss-pdb Viewer v4.0 (Arnold et al, 2006). 
1.4 Nucleotide sequence 

Two sequences of Bmp15 and Gdf9 gene from White 
goat in Guizhou province has been submitted to the 
GenBank database libraries and given the accession 
numbers of FJ429281 and FJ429282. 


2 Results 
2.1 Identification of Bmp15 and Gdf9 genes in White 
goat 


After comparing of DNA sequences encoding 
Bmp15 and Gdf9 genes, specific oligonucleotides were 


GTTTTACTTGACAGGAGTCTGTTAACGAC 


designed as primers to amplify fragments contained 
exons from genomic DNA of White goat. PCR products 
were cloned and sequenced (Fig. 1). Assembly sequences 
were obtained from repeated results of five prolific goats. 
Based on the blast by NCBI database, it contained 
complete exon! (328 bp) and exon2 (857 bp) in 
Bmp15-like gene, with partial 5'-flank region (215 bp) 
and intronl (38bp). Complete exon1 (397 bp) and exon2 
(965 bp) of Gdf9-like gene in White goat were screened 
as well, together with 5'-flank (50 bp), part of intron 1 
(23 bp) and 3'-flank region (21bp). They were 99% — 
100% identical to Bmp15 of Jining, Yunling and Boer 
goat breeds, and 99% to Gdf9 of Jining goat breed. The 
joined exon 1 and exon 2 of Bmp15- and Gdf9-like 
sequences encoded for 394aa of Bmp15- and 453 aa of 
Gdf9-like prepropeptide respectively by conceptual 
translation, which were very similar to that of other goat 
breeds. Highest peptide sequence identity was 98.7% — 
100% to that of Boer, Jining and Yunling goat breeds. It 
confirmed that the amplified fragments from White goat 
were BmpI5 and Gdf9 gene. The deduced precursor 
protein of BMP15 and GDF9 contained the proteolytic 
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Fig. 1 PCR products of exons fragments of Bmp15 and Gdf9 
in White goat 

Lane M: DL2000 DNA marker; Lane 1: exonl fragments of Bmp15 

(574 bp); Lane 2: exon 2 fragments of Bmp15 (862 bp); Lane 3: exon 1 

fragments of Gdf9 (461 bp); Lane 4: exon 2 fragments of Gdf9 (995 


bp). Lane 5: negative control. 
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processing site of RXXR in proregions to yield mature 
peptides, with 125 residues in BMP15 and 135 residues 
in GDF9. Each peptide contained a conserved 6-Cys 
domain which made up three pairs of disulfide bonds 
forming the cysteine-knots. In mature BMP15 peptide, 
seven amino acids changes were detected from the five 
female goats with D17G, PS50L, L65H, H70Y, V85A, 
P99I and G107R. It was interested that three amino acids, 
phenylalanine, serine or isoleucine, appeared at residue 
99 of the matured BMP15 (denoted as P99I/S99]). Fewer 
substitutions of amino acids were found out from the 
mature GDF9 peptide sequences, which were V18A, 
P78Q, V79I and T81A. Furthermore, GDF9 peptide 
sequence in mammals was much more conserved than 
that of BMP15. For example, the similarity of mature 
peptide between White goat and mouse was 87.4% in 
GDF9 while just 68.1% in BMP15. Taking the deposited 
sequences of BMP15 and GDF9 in White goat as 
represents, a multiple alignment containing amino acids 
sequences of different goat breeds was generated and a 
phylogenetic tree was constructed with other mammals 
and opossum, platypus, chicken and fish (Figs. 2, 3). The 
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Fig. 2 The phylogenetic tree of BMP15 prepropeptide by MEGA4.0 program using Neighbor-Joining method 


The bootstrap consensus inferred from 1 000 replicates. The percentage of replicate trees in which the associated taxa clustered together was shown 


next to the branches. The phylogenetic tree was linearized assuming equal evolutionary rates in all lineages. The tree is drawn to scale, with branch 


lengths in the same units as those of the evolutionary distances used to infer the phylogenetic tree. The evolutionary distances were computed using 


the Poisson correction method and Pairwise deletion. BMP15 prepropeptide (394aa) of White goat was much closer to the other three goat breeds 


(Yuling, Jining and Boer goats). The cited sequences in GenBank were as following, White goat: FJ429281; Jining goat: ABH08952; Yunling goat: 
ACF24868; Boer goat: ACF24873; sheep: AAF81688; cattle: AAS99651; pig: AAL58885; horse: XP_001496273; dog: XP_549005; human: 
AAC99768; mouse: AAC99766; rat: CAB41039; opossum: AAO12760; platypus: XP_001510645; chicken: AAU29611; zebrafish: AAI15086. 
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Fig. 3 Phylogenetic tree of GDF9 prepropeptide by program MEGA4.0 using Neighbor-Joining method 
All of correction and option was the same with that of BMP15. GDF9 prepropeptide of White goat was clustered together with other two goat 
breeds (Jining goat and goat described by Wu). The cited sequences in GenBank were as following, White goat: FJ429282; goat described by 
Wu: AAU09020; Jining goat: ABR10699; sheep: AAC28089; cattle: AAG38106; horse: XP_001504477; pig: AAT48523; human: AAH96230; 
mouse: AAH52667; rat: AAD 16406; opossum: XP_001372019; chicken: NP_996871; zebrafish: AAI08014. 


genetic relationships among sequences were assessed by 
MEGA4.0 program (Tamura et al, 2007) using 
method. BMPI5 and GDF9 
prepropeptide of White goat, Jining goat, Yunling goat 


Neighbor-Joining 


and Boer goat constituted a branch in the phylogenetic 
tree. 
2.2 Analysis of the divergence of BMP15 and 
GDF9 prepropeptide 

Higher identity of BMP15 prepropeptide sequence 
was presented between White goat and other mammals 
with percentage of 54.2%—99% while lower one was 
found between White goat and chicken (45.1%), platypus 
(40.7%) and fish (28.2%) (Fig. 2). In GDF9, the 
distances of White goat from other mammals were much 
closer than that from avian (58.3%) and fish (41.1%) 
(Fig. 3). 

Compared among vertebrates taking the sequence 
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of White goat deposited in the GenBank as template, the 
deletion mutations in BMP15 prepropeptide increased 
from mammals to fish (Fig. 4a, Fig. 5) and resulted in 
shorter BMP15 prepropeptide from mammals with 391 
—394 residues to chicken with 350 residues, in which 
deleted amino acids were 3 in human, 13 in opossum, 
and 45 in chicken. The BMP15 prepropeptide of 
zebrafish was reported as 394 residues which resulted 
from deletion of 19 amino acids and compensated by 
insertion of 18 residues. On the other hand, the length of 
GDF9 prepropeptide kept in 441—453 except for the 
length of zebrafish decreased to 418 residues, but both of 
deletion and insertion occurred more frequently from 
mammals, opossum, chicken to fish (Fig. 4b). 
2.3 The tertiary model of BMP15 and GDF9 in 
White goat 
The identity of mature peptide between BMP15 and 
b B Deletion [Ø Insertion 
20 


-10 
-20 
-30 
-40 
-50 


Deletion/insertion 
residues of GDF9 


Hum 
Mou 


Rat 


Goat 
Pos 
Pla 
Chi 
Fish 


Hor 
Dog 


She 
Pig 


Fig. 4 Profiles of deletion and insertion mutations in BMP15 (a) and GDF9 (b) prepropeptides 


The deletion numbers of BMP15 or GDF9 prepropeptide were showed downward column by negative values. Insertion mutation 


numbers produced in BMP15 or GDF9 prepropeptides standed upward by positive values. 













598 Zoological Research Vol. 30 
Fish RWHSCVRGAEHE 12 
White MVLLSILRILLLWGLVLFMEHRVOMT OVGOPSIAHL PEAPTLPLIQELLEEAPGKOQRKPRVL HGHPLRYMLELYOQRSAD 19 
Sheep cst veeses Nos cuesucssdvecase ses ereresescsisdserussesvtassasecnses Bo sacavronsudessse 78 
Cattle PROT ee TE Ee ey PLE IO Oh IC OCH ee TeMecccutccvesesses 79 
Pig overss if) SPCR Eee Te thee aCucaceRecesscactesess Q. Hac cccvscccccnces 79 
DOG. wmesstsrsosro T k E | cere N.LMADT.S Rie. ce cede ewes Q. He cccncceccvccces 18 
HOXGG Soran noe Kiveaes Riiseoss Sree LPA.Viseces Lis ove es ear E ASETETA 78 
HOMAM Sasccsrose es PE py ra BS eric Brives S. -S E T RTTEEEE Rissa 19 
Rat Beak A ot eee K...AK.EW. .TTL.A.H....SSLD.AK..... EM. QWPQH.Y...... | Ser 17 
Mouse eMacTessaee XaaVacest . -AKP.W..T.L.ADD....S.LD.AK..... EXM. QWPQH.Y...... K. HH... 1? 
Possum .GĦADTGxx.I..A. LFALRLSSPTGAHA. DAV. PAQVPTL.. -LGA..Q...RNSP.ROLTRH....Q...... ROW 7? 
Platypus A. .RE FLA. xxxxx, PW. LAFRGAGHPO . xxx, AR, SE. SEL... .GKT.xxxxx, SIFTH. Q. Q.. D.. £.. N 66 
Chicken -A..RPFTA. ..LTVL.SWAAxxxxxxxxxxxxxxxSQT.P...L. A. ‘RAQ. . .S. GURGGARS. Q... eseeoooo 65 
Fish GYORPHGLR. CVLSCLEVLHIF TRVAGHIAS PSHE GVAS TEVHRNRHPKRRN. HFRPPVESRTDDDGM. L.. S. RIA. 92 
White RE ge Sal RG STN A OLSEN” Ci 159 
SOG CO RiNeccad Sew SCREW AES RRMA CRU ERRHENUCCERE RRS COV A CRE RUA SSS RRS EMRE BY 158 
Cattle ETE E dcaia bale EOT oo NTA e SKUSE COR Clone CMOle Se SIE RMR RE RE ew AS v.: e 159 
Pig oRRa:6.0:0-2:6:6:016-5.06;0:0:6/0.0:0.0'0 WG. sses. P.sosoescescsososorsocacococoseosss DF.. secsoos Lees 159 
Dog esate wie wre'e-eere wien ove "ee MD ocie:tic:e i R ATEA ENT cee Essence b AT e s ccseseve I.Q. 158 
Horse M MASc ee KciWusscess 1 EEATT T P T N a P LA 158 
Human SU sie ste:a'e'0:5:6 nivevacarain's Sk Pere 1 Oe, rey Mere o AE A H. 0. BEB cease V..H 159 
Rat a AONA T E E A AT AS.E..... E RE ate e FRISII VP.C 157 
Mouse PH TE TET Reises K.S.HTV. Prinse Vicccons BS Osc Belesuessicveves VE PETET T.VP.C 157 
Possum QH.R...K.VER.NT...... AERR.L.P oc PoXV 2c wees acccs DDH... sA Ai PRBS. sccsesvcss Sa 157 
Platypus QE.L.W....L...S....PAS..TGIST..T.FV.S.EY.MSVVQERQI. . . PSV. RKSDTOFL .QMKSRIR.N 146 
Chicken HE.R. .RG.SLSWHT . QAASHGGO. WA.R. YV. P DR DME MELV. ‘A. POS. P. PRGR.L. A. .LLPAAK 145 
Fish -D.R.KQHKIF. SATI. ‘LOASTTEKHF . PxxxTSSD . QYTYTVKYELHD . LLDKL . KASFM. .RSPxxMSSRL.YICEA 167 
White PTHHF PSSGRGS PKPSLL PKTWIEMD IMEHVGOKLWHHKGRRVLRLREVCOGHHPRGSEVLEF WWHGTSSLDTVFLLLYF 237 
Sheep S HH 236 
Cattle 237 
Pig 237 
Dog 236 
Horse 236 
Human ace ee -Q. 236 
Rat Rik... SVS. A. R. N. THCIQ 234 
Mouse R.K.L...KS..S...PMS.A...1. . THCIQ ee ° oe « oll $ 235 
Possum xxxxxTILLG. GSPGEA. -ER. i: -LTNYIQ.QVOPQ...... HIQVR.. HIQERT . xIGLG. ROALAT. A..Y. 229 
Platypus EPHVREPEFEEASSQLYR.TA.Q. T.FTAYIR.R. QKAH. OFSVW. HWK. R. HHO. ARSSSRAP.R..AASEAP...... 224 
Chicken APAVLL . PTAP .RxxxxxxHG. A. A. . TPYLSPANSSSG. xxT.T. .HI . VRHHSGRAATAAP xxxxxP .PADP....FL 210 
Fish SVTSLONPLE .DRITMGPRSR. . `T. VTD. . ESXXXKD .HVSEFA. YW. TKPEHKR . VAHRKREPPOHH .RAPL , ` TFL 244 
White NDTOHSVOKTKPLPKGLKEF TEID PSLLLRRIR QARS TASEVESESREHDGEESINCSEHEEQVSEQUESIDHEELARHE 316 
Sheep RN 315 
Cattle ae 316 
Pig -RH.. . R. Ssa 316 
Dog -HK. . ° ° 315 
Horse -HK. L..R. os 316 
Human -HK. A.F..R. Rede ..T...DG. Sh. ‘TAS. SK.S. 314 
Rat D. .DHESAQA.L.AR.QE.L SPF .M.SV..TC....D..C..Q. 313 
Mouse -DHDRVOG.L.AR.QE.L.DRES.F.M.SV...C. ‘E. DASC..Q.... RT a 314 
Possum ‘HFK. . PRMELPELLVGDPEGA.LLS.2..V..V.PVR. .A. -Q.L. QT xxxxx ah able a ain eae Moc cbblas eases 304 
Platypus . HOGY .EVG. Pxxxxxxxxxxxxxx.H..V....xxxxKA. .S.L.LUKxxxxxF...R.LL.TSKS...TN.V...NK 281 
Chicken - RHXxxSGTLPEXxxxxxxxxxxxxxPR.SR.E..TLLUHDL. . x¥LRDA.GDKSD...RS.P.. A... see ene R 272 
Fish EENKHP .EWG. SF. PLSRPxxxxxx0RT, .SKKS.. .V.DI.HFKOGLSRVIK.N.K.YSSS .D .KD ce see's Via elk 316 
me SERGE EROS Fae LE EL VDONVPOPSCVEYKZVEISTLLTEANGSTEYKEYEGMIMISCTCR oa 

l Seid'cin'n ate 0.8: sielewhg.s 6.6: a:0 bhp uie/e wet b 0. 6b's'e'e's e106 66 b's ele Se deca sed t noes ule. 68a cect eM ier 
CAUEIO: Capano Eas M ins be eae nI 0 eles ale ee Sree iis vise cnt e ee Sas oes Cee en aie we aie ale init 394 
Fig A oe | Prrirerere rier err ret SENT RICE NURI CET REE ORE OEE ee Disvek 394 
Dog.. ee eieeves ENAT Ni tescav pees ni RAEADR PEE E caenes va p| ae tenn ne Diy een ea 393 
Hirse =. wlesccens r TYT Rissrcesssesos Siv.nelews'g’e PA AITE Dorsa cere ewe eeso Bivecncwas 394 
Haan see ccs nA ERTS POTTE TOET EsQeeocSs cRerssscuuccwe W Ms satis s0baue aes FERT 392 
Rat 8 aama E a AI ASA yee. A o ATE L E A A E 391 
MOUS  wasdirce Lili Bisrs enestis Se eve TS 66a xen ss NEL.M...... Dic siev IATL beieek 6 392 
Possum Sacowe T | eee te! E RSE CR vais. ciocc:c A Venice P E E Deseoso tao 382 
Platypus FN....Q.L...I.HI......... Mesura R.HRT..R.C Q@.R.. F..K.D..1V ERD. Vi osasse 359 
Chicken Rasin L.RD.YHA....VV bp aeMaweltsen's ee R.S... VM. GD, cccccove We ccMeccen 350 
Fish NLG..M.D...I.H. YY... 000s M.TFIS..GVADI.L........ KRM.M.VMGS..0.D..... nE EER 394 


Fig.5 Compare of BMP15 preproprotein in vertebrates with that of White goat as template 


The sign # denoted insertion, x as deletion and the dot (.) showed the same residue as template at the first line. The matured peptide of BMP15 


in White goat (125aa) was shaded in grey. All of the sequence of amino acid of BMP15 was the same as Fig. 2. 


GDF9 was 58.7%. In the comparison of human BMP7 
(hBMP7) with TGFB2, the three-dimensional structure is 
very similar even though their sequence identity is only 
36% (Griffith et al, 1996). Thus, in the light of crystal 
structure available for hBMP7 (pdb ID: cllx5A), the 
tertiary structures of mature BMP15 started from 22 to 
125 and GDF9 from 32 to 135 were built using Phyre 


program on line (Fig. 6, 7). In large part, the analogy to a 
curled left hand for TGFB2 (Daopin et al, 1992) was true 
in BMP15 and GDF9 monomers as well (Fig. 7a, 7b). 
The BMP15 or GDF9 cysteine-knots constituted the core 
of the monomer and consisted of three disulfide bonds: 
C24-C90, C53—C122 and C57—C124 in BMP15 and 
C34—C100, C63— C132 and C67—C134 in GDF9 (Fig. 
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6). Four strands of antiparallel B-sheet, which emanated 
from the knots, formed two fingerlike projections (Fig. 
7). Helix a2 (Fig. 6) located on the opposite end of the 
knots and lay vertically to the axis of two fingers thereby 
forming the heel of the left hand. The smaller finger 1 
contained Bl, al, B2 and $3, and the larger finger 2 
comprised of B4, B5, B6, and B7 (Fig. 6, 7). N-terminus 
of BMP15 and GDF9 corresponded to the thumb of the 
hand (not shown). The palm of the hand included the 


cysteine-knots and 31 residues contained in helix a2 (Fig. 


6, 7). The dimer of BMP15 and/or GDF9 was formed by 
a palm-to-palm orientation of the two monomers with an 
uncovalent binding (Fig. 7c). 


3 Discussion 


Bmp15 and Gdf9 genes are known to be specifically 
expressed in oocytes and to be essential for female 
fertility in sheep. Partial defects in Bmp15 and Gdf9 are 
associated with the increased ovulation rate in sheep 








(Galloway et al, 2000). We have investigated the 
naturally occurring mutation in Bmp15 and Gdf9 genes 
of White goat in Guizhou province and confirmed that 
the mutation of FecX® was present in White goat 
population (Lin et al, 2007). In addition, a different 
mutation from sheep was found in exon 2 of Gdf9 gene 
which was only detected from prolific White goat (Du et 
al, 2008). It was confirmed that heterozygous carriers 
contained one copy of mutation at 791 bp (G/A) in Gdf9 
exon 2 which resulted in the substitution of valine to 
isoleucine at residue 79 of GDF9 mature peptide (V791). 

To detect other diversities of BmpI5 and Gdf9 
genes, complete coding regions of exons in Bmp/5 and 
Gdf9 were cloned from five White goats of Guizhou 
province. Sequencing results showed that, like other 
TGFB superfamily members, Bmp/5 and Gdf9 of White 
goat were translated as prepropeptides composed of a 
signal peptide, a large proregion and a smaller mature 
peptide of 125 residues in BMP15 or 135 residues in 


gBMP15 1 QAGS | ASE VPGPSREHDGPES 21 Thumb 
gGDF9 1 DQESVSSELKKPLYPASA NLSEYFKQFLFPQ 31 
hRMP7 1 STGSKORSQNRS KTPKNGEAL RMANVAENSSSDOR 35 
fl al p3 
gBMP15 22 N QCSLHPFOVSFQGLGWDH WH LAPHLYTPNYCKGV C 57 Finger 1 
gGDF9 32 NECELHDFRLSFSGLKWON WI VAPHKYNPRY CKGD C 67 
hEMP7 36 GACKKH ELY VSFRDLGWQD WI LAPEGYAAYYCEGEC 71 
fil f2 al o2 f3 p4 p5 
o. 2 
gBMP15 58 PRVLYYGLNSPNHAI | QNLWNELVDONVP QP 88 Palm 
gGDF9 68  PRAVGHRYGSPVHTMYQN! | HEKLDSSWPRP 98 
hBMP7 72 AFPLNSYMNATNHAIVOQTLY HF | NPETYP KP 102 
a3 
f4 BS p6 P7 
gBMP15 89 SCVPYKYVPISILLIEANGSILYKEYEGMIAQSCTCR 125 Finger2 
gGDF9 99 SCVPAKYSPLSVLAIEPDGSIAYKEYEDMIATKCTCR 35 
hBMP7 103 CCAPTOLNAISVLYFDDSSNVILKKYRNMVVRACGCH 139 


p6 


p7 


Fig. 6 Structure-based sequence alignment of 6-Cys domains in White goat BMP15 (gBMP15), White goat GDF9 
(gGDF9) and human BMP7 (hBMP7, pdb ID: c1lx5A) 
Size of C (cysteine) in sequences was larger than the others in black font, while residues conserved in the mature peptides were shaded in 
grey. Sequences in B-sheets (81-87) were underlined and helixes (a1-a3) were dotted. The second structures in both gBMP15 and gGDF9 


(on the sequence of gBMP15) were similar while divergent from hBMP7 (under the sequence of hBMP7). Position of thumb, fingers and 


palm were denoted at the right side of the sequences. 
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Fig. 7 Homology-based model of BMP15 and GDF9 mature peptides in White goat 
Residue 22-125 of BMP15 (FJ429281) and 32-134 of GDF9 (FJ429282) were modeled on the basis of homology to the experimentally 
determined tertiary structure of human BMP7 (Griffith et al, 1996) with two fingers in blue and helix a2 in red located in the palm. 
Heterodimer model of BMP15 and GDF9 in White goat was constructed on the basis of homology to TGFB1 dimer (Greenwald et al, 2003). 


a: monomer model of BMP15 with substitutions at S32 (white balls) at finger 1, N66 (purple balls) at palm and $99 (green balls) at finger 2 


were highlighted. b: monomer model of GDF9 with substitutions at S77 (grey balls) and V79 (orange balls) were distinguished. c: dimer of 
gBMP15 (up) and gGDF9 (down). Diversity residues of BMP15 were highlighted as: P50 (could not be seen from the outside of the dimer 


as it located inside of finger 1), L65 (green balls, inside of the palm), H70 (dark green balls, surface of the palm), V85 (yellow balls, near to 
H70), S99 (white balls), G107 (orange balls). Polymorphism sites of GDF9 were showed as: S77 (yellow balls), P78 (orange balls), V79 


(green balls) and T81 (white balls, inside of the palm). The position of receptors was denoted according to McNatty et al (2004). 


GDF9. Deduced mature peptides of BMP15 and GDF9 
contained cysteine-knots which were conserved in TGFB 
superfamily. Apart from the detected polymorphisms in 
matured BMP15 (S991) and GDF9 (V791) (Du et al,2008; 
Lin et al, 2007), seven and three substitutions were found 
in mature peptide sequences of BMPI5 and GDF9 
respectively in White goat population. And another two 
changes of BMP15 (S32G and N66H) were found from 
that of Jining and Yunlin goat breeds. The effect of these 
polymorphisms in BMP15 and GDF9 on their function 
would be an interesting issue to address in the future and 
some clues from their structure analysis might be helpful. 

The structure of BMP15 or GDF9 peptide has not 
yet been determined by experiment assay. Based on 
human BMP7 (Griffith et al, 1996), 
three-dimension of BMPI5 and GDF9 
monomers of White goat were modeled by computer on 


structure 
structures 


line (Fig. 7). Each of them could be compared to a left 
hand with a palm and two fingers. To date, it still 
remains to be elucidated that these molecules exist as 
monomers or dimers in vivo. Evidences from TGFB 
members show that for the proteolytic cleavage of 
prepropeptides to occur, the prepropeptide must first 
dimerize (Hogan, 1996). BMP15 and GDF9 lack a 
cysteine which is conserved in TGFB superfamily 
disulphide 
(McPherron et al, 1993). However, in vitro studies using 


required for the  intersubunit bridge 


cell expression systems indicate that both BMP15 and 


GDF9 can form non-covalent homodimers and 
BMP15/GDF9 heterodimers (Liao et al, 2003). Therefore, 
a heterodimeric BMP15 and GDF9 of White goat was 
modeled in the light of TGFB1 dimer (Hinck et al, 1996). 
GDF9 and BMP15 molecules linked non-covalently by 
palm contacted each other, and four receptors binding 
around the dimer (Fig. 7c). 

In these seven polymorphism sites of BMPI5 in 
White goat individuals, D17G was on the thumb of the 
molecule which could not be proposed as the thumb 
structure needed further experimental data. PSOL located 
inner part of finger! and L65H was inside of the palm. 
Both of them might not affect the stability of dimer 
and/or the affinity of the dimer binding to its receptors. 
H70Y was at the side part of the palm and it might have 
little role on the function of BMP15 for it was far away 
from the connection region of the dimer and its receptors 
(Fig. 7c). V85A was at the connection region of the 
dimer which might be unimportant because both of 
valine and alanine are hydrophobic. G107R was at the tip 
of finger 2 which was the connection region of dimer 
with receptor type 2 like BMPR-II. S99I/P99I was 
between the seam of finger 1 and finger 2 and protruded 
from the back of the dimer which was important for the 
dimer binding to its receptor 2 (Fig. 7c). This coincides 
with previous reports that S99I in sheep (FecX”) affects 
receptor BMPR-II binding by structural determination of 
ligand-receptor complexes with BMP2, BMP7, and 
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activin A, and by site-directed mutagenesis (Greenwald 
et al, 2003; Thompson et al, 2003). For the difference of 
BMP15 between White goat in Guizhou province and 
Jining, Yunlin goat breeds, S32G and N66H were found 
at the surface of fingerl and palm of BMP15, 
respectively (Fig. 7a). The ovulation of sheep is potently 
inhibited when immunized by directed against the 
peptide fragments contained S32 and N66 residues of 
BMP15 (McNatty et al, 2007). Thus, like S99, S32 and 
N66 might participate in the binding of BMP15 with 
receptors but not the formation of dimer. 

For GDF9 of White goat, like D17G in BMP1S, 
V18A was in the thumb of the handed-like molecule, 
which was also observed in goat Gdf9 described by Wu 
et al. (GenBank no. AAU09020). T81A was hidden in 
the palm of GDF9 (Fig. 7c) which might be unimportant. 
However, two sites, P78Q and V79I, were adjacent to 
S77 which was FecG" mutation (S77F) found from 
sheep (Fig. 7b). All of these three residues occupied the 
side part of the GDF9 palm. Previous finding has showed 
that FecG” mutation changes an uncharged polar serine 
residue to a nonpolar phenylalanine in a region of the 
binding to ALKS5 (type I receptor of TGFB superfamily) 
(Kirsch et al, 2000). This change at S77 also affects the 
H80 in helix a2 of the mature GDF9 peptide. In BMP7, 
this conserved histidine H80 is important in dimmer 
binding which exhibits hydrogen bonding to three 
residues of the homodimer of BMP7 (Griffith et al, 1996) 
and heterodimer of TGFB3 (Mittl et al, 1996). On the 
other hand, P78Q changed from nonpolar proline to polar 
glutamine, and the diversity at residue 79 from valine to 
isoleucine increased a mehyl group at the side chain of 
amino acid. Further, site I79 was detected only from the 
prolific population of White goat (Du et al, 2008). It was 
reasonable that polymorphism sites, P78Q and V79I, 
might have an influence to the GDF9 binding to 
receptor! (Fig. 7c) such as ALKS. 

Deduced amino acids sequences of BMP15 and 
GDF9 in White goat were further compared with that of 
other animals. The length polymorphism presented in 
both of BMP15 and GDF9 prepropeptides. Compared 
with GDF9, higher frequency of deletion mutations was 
found inthe prepropeptides of BMP15 in opossum, 
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